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semiclassical  limit  of  the  nonlinear 
schrodinger-poisson  equation 
with  subcritical  initial  data  * 

HAILIANG  LlUt  AND  EITAN  TADMORt 

Abstract.  We  study  the  semi-classical  limit  of  the  nonlinear  Schrodinger-Poisson  (NLSP)  equa- 
tion  for  initial  data  of  the  WKB  type.  The  semi-classical  limit  in  this  case  is  realized  in  terms  of 
a  density- velocity  pair  governed  by  the  Euler-Poisson  equations.  Recently  we  have  shown  in  [ELT, 
Indiana  Univ.  Math.  J.,  50  (2001),  109-157],  that  the  isotropic  Euler-Poisson  equations  admit 
a  critical  threshold  phenomena,  where  initial  data  in  the  sub-critical  regime  give  rise  to  globally 
smooth  solutions.  Consequently,  we  justify  the  semi-classical  limit  for  sub-critical  NLSP  initial  data 
and  confirm  the  validity  of  the  WKB  method. 


1.  Introduction.  Consider  the  nonlinear  Schrodinger-Poisson  (NLSP)  equation, 

=  -  yAf  -  |2  -  c)]f ,  x  G  IRd,  fG  IR+,  (1.1) 

subject  to  the  initial  condition 

ipo(x)  =  V/ftWeiSo(x)/£,  *  G  IRd.  (1.2) 

Here,  ^>e(-,t)  is  a  complex- valued  wavefunction  depending  on  the  scaled  Planck  con¬ 
stant  e  ~  H,  with  k  and  c  >  0  being  scaled  physical  constants.  The  NLSP  equation 
(1.1)  has  been  studied  in  different  contexts,  and  in  particular,  as  the  fundamental 
equation  in  semiconductors  applications,  with  c  >  0  standing  for  a  constant  back¬ 
ground  charge  and  k  ~  A~2,  A  <<  1  being  the  Debye  number.  We  refer  the  reader  to 
the  recent  review  [GLM]  and  references  therein. 

The  connection  between  Schrodinger  equations  and  the  classical  hydrodynamics 
was  already  noted  in  1927  by  Madelung,  in  the  context  of  semi-classical  limit  of  the 
Nonlinear  Schrodinger  (NLS)  equation.  To  this  end,  one  identifies  two  physically 
relevant  observable  quantities  —  the  fluid  density  pe  :=  \ipe\2,  and  the  fluid  velocity 
ue  :=  e\7 xargil>e .  Indeed,  introducing  Se  as  the  phase  of  the  wave  function 

=  \J pe(x,  t)eiS‘(x,t^e 

and  separating  real  and  imaginary  parts  in  the  NLSP  equation  (1.1),  one  obtains  the 
following  irrotational  flow  equations  for  the  density-velocity  pair,  ( pe,ue  :=  \7xSe), 
consult  e.g.,  [LL] 


p\  +  v  •  (p£«e)  =  0) 


u\  +  ue  ■  Vtte  =  k\7Ax  1(pe 


'A# 

.  Vf  _  ■ 


(1.3) 

(1.4) 


The  resulting  system  amounts  to  a  dispersive  regularization  with  the  square  of  the 
Planck  constant  playing  the  role  of  the  amplitude  of  dispersion.  If  we  argue  formally 
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that  this  0(e2)  term  on  the  right  of  (1.4)  is  negligible  as  e  — >  0,  then  the  corresponding 
limiting  system  becomes  the  Euler-Poisson  equations 

Pt  +  V  •  (pu)  =  0,  (1.5) 

ut  +  u  •  Vu  =  fcVA“1(p  —  c),  (1.6) 


subject  to  given  initial  data 

p(x,  0)  =  po(x),  u(x,0)  =  VS0(x). 

The  dispersive  regularization  term  in  (1.4)  represents  a  quantum  correction,  and  the 
corresponding  system  (1.3)-(1.4)  governing  the  observables  ( p€ ,  ue)  could  be  viewed  as 
Quantum  Euler-Poisson  (QEP)  system,  a  specific  example  to  larger  class  of  quantum 
hydrodynamical  systems,  [SMM]. 

The  above  argument  is,  of  course,  only  formal.  The  Madelung’s  transformation 
relies  on  the  assumption  that  the  amplitude  pe  does  not  vanish  and  that  the  phase 
Se  remains  nonsingular,  for  otherwise  the  transformation  is  not  well-defined  and  the 
equations  (1.3),  (1.4)  become  singular,  even  though  NLSP  itself  is  still  regular.  The 
challenging  issue  here  is  to  justify  such  a  dispersive  limit.  The  asymptotics  of  ob¬ 
servable  quantities  as  e  — >  0  is  known  as  ‘semi-classical’,  expressing  the  passage  from 
quantum  to  Newton  mechanics,  where  the  time  and  distance  scales  become  large 
enough  relative  to  Planck’s  constant.  The  ‘semi-classical’  limit  problem  associated 
with  the  NLSP  equation  (1.1)  can  be  summarized  as  follows.  We  consider  the  family 
of  solutions  for  the  NLSP  Cauchy  problem  (1.1),  {tpe,  e  >  0},  subject  to  the  WKB-type 
initial  data 


tpe{x,0)  =  V/M£)eiSo(x)/€,  x  G  IRd  (1.7) 

where  po  and  Sq  are  given  smooth  functions.  The  purpose  is  to  determine  the  behavior 
of  the  wavefunction  if)*  as  e  — >  0,  through  the  dynamics  governing  semi-classical  limits 
of  its  observables,  the  density,  p  =  lim  |?/A|2,  and  the  velocity  u  =  lim  eVa,arg?/;e. 

As  noted  earlier,  this  argument  is  self-consistent  as  long  as  the  solution  of  the 
Euler-Poisson  system  (1.5)-(1.6)  remains  classical.  Motivated  by  the  recent  regularity 
results  for  general  a  class  of  isotropic  Euler-Poisson  equations,  [ELT],  we  turn  to 
justify  the  (strong)  semi-classical  limit  of  the  NLSP.  In  this  paper  we  restrict  our 
attention  to  the  one-dimensional  problem 

=  -yV’L  +  VL  =  KW\2  -  c).  (1.8) 

The  corresponding  one-dimensional  QEP  system  reads 
Pet  +  ( Peue)x  =  0, 

v*  +  u*u%  =  VZ+^(y^}  ,  V*x  =  k(p€  —  c), 

and  we  are  inquiring  whether  the  one-dimensional  NLSP  (1.8)  has  a  ‘faithful’  de¬ 
scription  in  terms  of  its  formal  semi-classical  limit  system,  the  one-dimensional  Euler- 
Poisson  equations, 


(1.9) 

(1.10) 


pt  +  ( pu)x  =  0, 

ut  +  uux  =  Vx,  Vxx  =  k(p  -  c) . 


(1.11) 

(1.12) 
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In  [ELT]  we  have  shown  that  the  one-dimensional  Euler-Poisson  system  (1.11)- 
(1.12)  with  initial  data,  (po,  uo),  f  (po(’)  —  c)  =  0,  below  certain  critical  threshold, 
i.e., 

\dxu0\  <  k(2po  -  c),  (1.13) 

retains  a  globally  smooth  solution.  Equipped  with  this  regularity  result  for  the  limit¬ 
ing  Euler-Poisson  system,  it  is  natural  to  inquire  about  the  regularity  of  its  dispersive 
regularization,  the  QEP  equations  (1.9)-(1.10).  Our  main  results  show  the  existence 
of  globally  smooth  solution,  ip€,  of  (1.8),  which  assume  sub-critical  smooth  initial  data, 

(1-7), 

\d2xS0\  <  v/fc(2|^0|2-c). 

We  then  justify  the  semi-classical  limit  for  the  observables  (pe,ue)  connected  with 
these  global  solutions  and  we  conclude  by  verifying  their  WKB  expansions. 

We  comment  in  passing  on  two  related  difficulties  associated  with  justifying  the 
global  semi-classical  limit  in  this  case.  First,  the  QEP  equations  represent  a  dispersive 
regularization  of  the  limiting  Euler-Poisson  equations,  and  as  such,  they  do  not  lend 
themselves  to  the  standard  energy  method.  Instead,  the  analysis  of  small  dispersive 
limits  usually  proceeds  by  exploiting  sufficiently  many  integrals  of  the  motion.  A  case 
in  point  is  the  recent  study  of  semi-classical  limit  of  the  defocusing  cubic  NLS  carried 
out  in  [JLM2] .  In  our  case,  however,  the  nonlinearity  of  the  dispersive  correction  yields 
only  two  such  integrals  that  we  know  of  —  the  total  density  and  energy  recorded  in 
(2. l)-(2.2)  below.  Second,  the  regularity  of  the  Euler  Poisson  equations  with  sub- 
critical  initial  data  is  based  on  tracing  along  particle  paths,  the  delicate  balance 
between  the  focusing  effect  of  convective  terms  and  the  repulsive  Poisson  forcing.  We 
lack,  however,  the  usual  energy  dissipation  estimates  in  this  case.  Consequently,  one 
should  not  expect  to  extend  a  classical  energy  approach,  whose  general  framework  in 
the  context  of  quantum  hydrodynamics  can  be  found  in  [JLM1],  in  order  to  justify 
the  semi-classical  limit  in  the  present  setting  of  sub-critical  NLSP  problem. 

To  circumvent  these  difficulties,  our  analysis  of  the  2x2  QEP  system  (1.9)-(1.10) 
is  carried  out  in  terms  of  a  modified  Madelung’s  transformation  of  the  form,  consult 
[Gre], 


i>e(x,t)  =  Ae(x,t)el*^x’t)/e. 


Here  Ae  is  a  complex  valued  amplitude  which  takes  into  account  the  slowly  varying  part 
of  the  solution  and  represents  the  fast  scale.  There  is  a  one-to-one  correspondence 
between  these  modified  variables,  (qe  :=  $ ex,Ae  :=  ae  +  ibe)  and  the  observables 
( pe,ue  =  S%.)  connected  with  the  usual  representation,  ipe  =  yfpf exp{iSe / e} , 


Pe  = 


ue  =  qe  +  e 


aeb%  -  bea[ 

2  V(o£)2 


:=  Si 


=:  T! 


(1.14) 


This  leads  to  an  equivalent  3x3  modified  system,  (2.8)-(2.10),  governing  the  modified 
velocity  gradient  qe  and  a  slowly  varying  complex  valued  amplitude  Ae  =  ae  +  ibe. 
Unlike  the  usual  WKB  variables,  the  governing  equations  for  Ue  :=  ae,  6e]  form  a 
strongly  hyperbolic  system,  which  enables  us  to  bridge  between  the  associated  NLSP 
solution,  Ue ,  and  its  semi-classical  limit  solution  of  the  Euler-Poisson  equations. 
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The  rest  of  the  paper  is  summarized  into  main  four  theorems.  We  begin  in  §2  with 
the  following  local  existence  statement. 

Theorem  1.1  (Local  existence).  Consider  the  NLSP  equation  (1.8)  subject  to 
a  sequence  of  WKB  initial  data,  ipo(x)  =  A()(x)exp{i'hQ(x) /e},  with  (<3>g(a;),  Ag(a;)) 
uniformly  bounded  in  (7J4(IR), H2(\R)).  Then  there  exists  To  >  0, 

T0  ~  l/(||$ol|ff4(IR)  +  P5||*2(|R)), 

such  that  (1.8)  admits  a  solution  of  the  form 

ipe(x,t)  =  Ae(x,t)exp{i$e(x,t)/e},  0  <  t  <  T0. 

Moreover  (Ae(x,t),$£(x,t))  remain  bounded  in  T°°([0,  T0],H2  (IR))  x 
T°°([0,To],iL4(IR))  uniformly  in  e. 


Theorem  1.1  tells  us  that  the  modified  variables,  U£  :=  [d£  :=  $£x ,  a£,  b£]  converge 
to  a  limit,  U  :=  [d,  a,  b\.  This  limit  is  governed  by  a  modified  limit  system,  (3.1)-(3.3), 
which  is  studied  in  §3.  Our  next  result,  outlined  in  §4,  complements  Theorem  1.1 
with  an  error  estimate  for  \\U£  —  U\\. 

Theorem  1.2  (Error  estimate).  Consider  the  NLSP  equation  (1.8)  subject  to 
a  sequence  of  WKB  initial  data  ipoix)  =  x4g(cc)ea7^*{zcI>Q (ar) / with  (4>o(;r),  Aq(x)) 
uniformly  bounded  in  (TT4 ( IR) ,  TT2 (IR) ) .  Let  U£  =  [d£  =  <&%x,  a£,  b£],  t  <  To  denote 
the  modified  variables  associated  with  the  local  NLSP  solution  A£(x,t)exp{i( I»£/e},  t  < 
To,  and  let  U  =  [ d,a,b }  denote  the  solution  of  the  modified  limit  system,  consult  (3.1)- 
(3.3)  below.  Then  H2  —  lim  U£(-,  t)  =  U(-,t),  t  <  T0,  and  the  following  error  estimate 
holds 


\Ue(-,t)  -  U(-,t)\\H2  <  Const.  \\Uo(x)  -  U0(x)\\H2 


t  <  Tg. 


The  modified  limit  system  (3.1)-(3.3)  is  equivalent  with  Euler-Poisson  equations 
(1.11)-(1.12).  In  particular,  the  IL2-regularity  of  its  limit  solution,  U(-,t )  is  dictated 
by  the  regularity  of  the  Euler-Poisson  equation.  We  now  invoke  the  global  regularity  of 
Euler-Poisson  solutions,  established  in  [ELT]  under  the  sub-critical  condition  (1.13). 
It  follows  that  ||C7(-,  t)||jj2  and  hence  \\Ue(-,  f)||#2  remain  bounded  for  t  <  Const.ln\ 
and  we  can  therefore  extend  the  local  solution  constructed  above.  The  argument  is 
classical  -  consult  e.g.,  [Co]  in  the  context  of  global  Navier-Stokes  solutions  during  the 
regularity  interval  of  Euler  equations.  We  conclude,  in  §5  with  the  following  statement 
of  global  existence. 

Theorem  1.3  (Global  existence).  Consider  the  NLSP  equation  (1.8)  subject  to 
WKB  initial  data,  ^(x,Q)  =  Ao(x)exp{i$o(x) /e} >  with  4>g(a;)  =  S'o(x)  €  TT4(IR),  and 
A0(x)  =  \/po  ( x )  €  H2(\R),  f  (po(-)  —  c)  =  0.  Assume  that  the  initial  data  satisfy  the 
sub-critical  condition 


|92S'o(o;)|  <  yjk(2p0(a)  -  c),  a  £  IR. 

Then  for  any  T  >  0,  there  exists  eo  =  £o (T)  such  that,  for  every  0  <  e  <  eo,  the  NLSP 
equation  (1.8)  admits  a  global  solution  of  the  form  =  A£(x,t)exp{iA>e(x,t) / e\  for 
all  t  <T,  where  $>l,x,Ae  =  ae  +  ib£  are  bounded  in  T°°([0,  T],  TT2 ( IR) )  uniformly  in  e. 
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Moreover,  its  semi-classical  limit  ( d  =  &xx,a,b)  is  the  global  solution  of  the  modified 
system  (3.1)-(3.3)  and  the  following  error  estimate  holds 

||^*M)  -  \\m  +  \\Ae(-,t)  -  ||jj2  <  ConstT-e,  t  <  T.  (1.15) 


Finally,  we  are  ready  to  convert  the  last  global  statement,  yielding  a  complete 
global  characterization  of  the  semi-classical  limit  of  the  one-dimensional  NLSP  in 
terms  of  the  usual  observable  variables,  pe,ue.  Relations  (1.14)  and  Theorem  1.3 
imply  the  H2  convergence  of  pe  — >  p{x,t)  of  order  0(e).  Moreover,  we  can  use  the 
fact  that  i71(IR)  is  an  algebra  to  conclude  that 


aebx  —  beax 

2iyV)2  +  Wf 


<  Const.e 

H1 


and  hence  the  Hi  convergence  of  ue  =  qe  +  0(e )  — >  u  of  order  0(e)  follows.  We 
conclude 

COROLLARY  1.1  (Global  convergence  revisited).  Consider  the  NLSP  equa¬ 
tion  (1.8)  subject  to  WKB  initial  data  (1. 7),  ifo(x)  =  \J po(x)exp{iSo(x) /e },  with 
(po(x),  So(x))  bounded  in  (i72(IR),  i74(IR)),  f(po(m)  —  c)  =  0.  Assume  that  the  initial 
data  satisfy  the  sub-critical  condition 


|S2S'o(Q!)|  <  \J k(2p0(a)  -  c),  a  G  IR. 


Then ,  for  e  small  enough,  (1.8)  admits  a  global  solution  of  the  form 


=  \J pe(x ,  t)exp{iSe(x,  t)/e} 

for  all  t  <  Const.ln-,  with  the  following  semi-classical  limit,  ( pe,  ue  :=  Sx )  — >  (p,  u) 
where  ( p,u )  is  the  solution  of  the  Euler-Poisson  equations  (1.5)  (1.6).  Moreover,  the 
following  error  estimate  holds 

II Pe{‘,t)  ~  p(-,t)\\H2  +  || ue(-,t)  -  u(-,t)\\Hl  <  ConstT  ■  e,  t  <  T  ~  In 


The  above  results  enable  us  to  justify  the  WKB  expansion. 

Theorem  1.4  (WKB  expansion).  Consider  the  NLSP  equation  (1.8)  subject  to 
sub-critical  WKB  initial  data,  if(x,  0)  =  A(x,Q)exp{iA>{x,Q) / e\ 

|^d>(a,0)|  <  \A(2|7L(a,0)|2-c). 

Fix  an  integer  N  >  0,  and  assume  that  [<f>a,a:(-,  0),  A(-,  0)]  G  Hs ,  s  >  2 N  +  2.  Then 
(1.8)  admits  a  solution  of  the  form  =  Ae exp{i<S>£ / ^e} , 

N 

A\x,t ))  =  J2Aj{x,t)C  +o(eN),  Aj  G  Hs~2j{\R), 

3=  0 
N 

=  +  o(eN), 

3= 0 


g  Hs~21+2( IR). 
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This  then  yields  the  usual  N -terms  WKB  expansion 

N  _ 

Ae(x,  t)exp{i$e(x,  t)/e}  =  ^  \J pj (x,  t)e^ exp{iSj (x,  t) / e}  +  o(eN),  t  <T  ~  Nln- 

j=o  e 

with  Sj  and  pj  determined  by  the  WKB  method. 

We  conclude  our  Introduction  by  mapping  yet  another  possible  route  to  justify  the 
semi-classical  limit  in  the  current  setting  of  NLSP.  In  the  first  step,  one  justifies  the 
passage  from  one-dimensional  NLSP  to  the  limiting  one-dimensional  Vlasov-Poisson 
equation, 

dtf(x,  £,  t)  +  £dxf( x,  £,  t)  +  Vxd^f(x,  £,  t)  =  0;  Vxx(x, t)  =  k(J  f(x,  £,  t)d£  -  c) . 

This  passage,  administered  through  the  Wigner  transform,  fe  :=  (Wi^e)(x,  £,  t),  was 
rigorously  justified  in  [LP]  and  [MM]  for  multidimensional  problems  in  the  mixed  type 
case  and  was  recently  proved  in  the  one  space  dimension  for  the  general  case  in  [ZZM]. 
The  second  step  requires  a  rigorous  justification  of  the  passage  from  Vlasov-Poisson  to 
Euler-Poisson  equations.  This  passage  is  administered  through  the  measure  solution, 

/( x,  £,  t)  =  p(x,  t)6(£  -  u(x,  t)), 

whose  first  two  moments  form  a  renormalized  solution  of  the  Euler-Poisson,  consult 
[GM],  In  particular,  for  sub-critical  (p0,  u 0),  the  Euler-Poisson  solution  ( p(-,t),u( •,  t)) 
remains  smooth,  [ELT]  and  it  coincides  with  such  renormalized  solution,  [GLM,  Corol¬ 
lary  3.2]. 

2.  Local  convergence  in  time.  We  can  show  that  there  are  two  global  invari¬ 
ants  for  the  quantum  Euler-Poisson  equation 

M(t)  =  [  (pe-  c)dx  =  M( 0),  (2.1) 

J IR 

m  :=  £  Qpe(tie)2  +  \{V*Y  +  e2((v^)x)2)  dx  =  E{ 0).  (2.2) 

Following  Grenier  [Gre] ,  we  use  a  modified  Madelung’s  transform 

^  =  Ae{x,  t)ei4>e(x,t)/e 

with  a  complex-valued  ‘amplitude’,  Ae,  representing  the  slowly  varying  part  of  the 
solution  and  an  argument  <I>e  dictating  the  fast  part  of  the  solution.  In  the  present 
context,  this  modified  transformation  plays  a  key  role  in  the  derivation  of  global  (in 
time)  hyperbolicity  which  in  turn  justifies  the  global  semi-classical  limit.  The  relation 
with  the  usual  observables  is  revealed  upon  the  change  of  variables 

=  y/feiSe/e  =  \Ae\exp{i(Se  +  earg  Ae)/e} 

which  yields 

:=  l^l2  =  |Af ,  (2.3) 

u*-.=  S*=<f  +  e{A*A%-A*%)l{2ip*),  qe  :=  K-  (2-4) 
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Substitution  of  ipe  =  Ae exp{i$£ / e}  into  the  NLSP  (1.8)  leads  to 


7  F  1  9 

-ieA\  +  ^  -  -A%x  -  ie^xA£  -  -A^%x  +  -A£(<1>£)2  -  =  0 


with 


Vxex  =  |^|2  -  c. 

We  split  the  above  equation  between  the  slowly  varying  part  governed  by  the  complex 
valued  A£,  and  the  fast  scale  governed  by  <1>C, 

*et  +  l(K)2  =  Ve,  V'x  =  k(\Ae\a-c), 

1  If 

At  +  *%A%  +  -A**lx- -A£xx  =  0. 

Setting  the  velocity  q£  =  4>£  and  the  electric  charge,  ee  :=  Vx/k,  gives 

qt  +  q£q£x  =  ke£,  e£  =  |A£|2  -  c,  (2.5) 

1  if 

Al  +  cfA%+-A£q%=-Alx. 

Separating  the  real  and  the  imaginary  part  of  Ae  :=  ae  +  ibe ,  the  amplitude 
equation  reads 


<  +  <fa%  +  \a£qex  =  ~^bexx,  (2.6) 

bt  +  <fb%  +  =  2a%x-  (2-7) 

To  justify  the  uniform  estimates  of  4>e  and  Ae,  the  velocity  is  dealt  in  terms  of  its 
derivative,  de  :=  qx  =  $%x,  satisfying  the  transport  equation 

dl  +  qedx  +  (cf )“  =  kex  =  fc(|ae  +  ibe\2  —  c). 

We  arrive  at  the  following  system 

det  +  qedx  =  k(\ae  +  ibe\2  -  c)  -  (de)2,  (2.8) 

at  +  qeaex  = -^aede  -  ^bexx,  (2.9) 

bl  +  qebx  =  -\bed£  +  (-a%x.  (2.10) 

We  observe  that  even  when  the  modified  variables  assume  initial  data  with  a 
real  amplitude,  (ae(x,  0),  be(x,  0))  =  0),  the  evolution  of  (2.10)  still  develops  an 

imaginary  component,  be(-,t )  ^  0. 

The  system  (2.8)-(2.10)  can  be  written  in  the  compact  form  in  terms  of  the  3- 
vector  Ue  :=  [de,  ae,  fr£]T 


e 

2 


JU1 


Q(U£,U£),  J  := 


0  0  0 

0  0-1 

0  10 


Ut  +  qeUex 


(2.11) 


SCHRODINGER-POISSON  EQUATION  WITH  SUBCRITICAL  INITIAL  DATA  477 


with  a  quadratic  nonlinear  term  given  by 


Q{Ue,Ue) 


k(\ae  +  ibe\2  —  c)  —  ( de )2 
-aede/ 2 
-bede/2 


Observing  that  (2.11)  is  a  symmetric  hyperbolic  system,  will  enable  us  to  proceed 
in  the  usual  fashion  with  classical  energy  estimate.  We  use  (-,  •),  ||  •  ||  to  denote  the 
usual  L2  products  and  norms,  and  we  seek  a  short  time  uniform  bound  of  ||[/£||^  := 
El<m\\dlxU%2.  Differentiation  of  (2.11)  with  respect  to  x  leads  to 


+  d™{<fu*x)  =  -J(d™U%x  +  d?Q{U*,  t/£),  m  =  0,1, 2. 


Taking  the  inner  product  of  the  above  equation  with  d™Ue  with  m  =  2  yields 

\jtWdlueW2  =  (dtdlu\dlu*) 

=  (  -  d2xw)  +  e-  (J(d2Ue)xx,  d lu<)  +  (d2Q(U\  C/£),  d2U £) 


=:/i 


For  the  first  term  we  have 


h  =  (-  d2x(q*U%,d2xU<)  =-t(l)  {dkxqedl~kU\d2xW)  =:  £  Ilk 

b— n  \  /  b—n 


2 


k—0 


k—0 


Integration  by  parts  yields  |/i0|  =  |  (qedxUe,  d2Ue)  \  =  \  \  (qxdxU€,  d2Ue)  \  < 
||d£||oo||C/e|||.  For  the  other  two  terms  we  have 

|/11|  =  2|(^C/£,^[7£)|<|K||00||C/£||i 

l^ial  =  \(qexxdxUe,d2Ue)\  <  ||4lloo||t/e||!, 

and  together  with  one-dimensional  Sobolev  estimate,  the  last  three  upper  bounds 
imply  that 

| Jr |  <  Const.(\\d€ \\oo  +  || ^ ||oo ) II ||  1  <  Const.\\Ue\\l  (2-12) 


Integration  by  parts  and  the  anti-symmetry  of  J  show  that  the  dispersive  term  van¬ 
ishes, 

h  =  ^(J(d2xu%x,d2u £)  =  -e-(Jd3xu*,d3xu<)  =  o. 

We  are  left  with  the  last  quadratic  term 

h  =  {d2xQ(u*,uo,d2xw) 

=  {{k{\a*+iV\2-c)-{d02)xx,d%x)  \{{a^)xx,a%x)  \{{b^)xx,b%x) 

=  '■  J3I  +  J32  +  J33- 


Some  further  calculations  imply  that  each  of  the  last  three  expressions  is  upper 
bounded  by  Const. ||  || i;  f°r  example 


'31 


<  ^ 


\K 


>11  dex 


\a%x\\  <  Const. \\Ue 
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with  similar  estimates  for  /32  and  J3 3.  The  estimate  for  ||<9™£/m||  with  m  =  0,1  can 
be  similarly  obtained.  A  combination  of  the  above  facts  yields  the  usual  Ricatti’s 
equation 


jtm\2<C0\\U*\\i  (2.13) 

The  energy  bound  independent  of  e  asserted  in  Theorem  1.1  follows  for  a  local  time 
interval,  [0,T0], 

3.  The  modified  limit  system.  To  study  the  global  convergence  rate  of  the 
QEP  system  (1.9)-(1.10),  we  consider  the  modified  Madelung  transformation  which 
led  to  the  system  of  three  equations  (2.8)-(2.10).  The  corresponding  3x3  limit  system 
takes  of  the  form 

dt  +  qdx  =  fc(|a  +  ib\2  —  c)  —  d2, 
at  +  qax  =  ad , 

bt  +  qbx  =  bd , 

with  q(-,t)  standing  for  the  primitive  of  d(-,t)  satisfying 

qt  +  qqx  =  ke,  ex  =  \a  +  ib\2  -  c.  (3.4) 


(3.1) 

(3.2) 

(3.3) 


This  modified  system  shares  the  same  regularity  properties  as  the  2x2  limit  Euler- 
Poisson  system,  (1.11-(1.12).  This  is  the  content  of 

Lemma  3.1  (Regularity).  Let  T*  be  the  maximal  time  interval  of  regularity  for 
the  solution  U  =  ( d,a,b)T  of  the  modified  system  (3.1)-(3.f).  Let  T  be  the  maximal 
time  of  regularity  for  the  solution  (p,  u)T  of  the  Euler-Poisson  equation  (1.11)-(1.12), 
subject  to  initial  conditions  (po  =  a2  +  &q,uo)t.  Then  we  have  T  =  T* . 

Proof.  Assume  (d,  a,  b)T  is  a  smooth  solution  of  the  modified  system  (3.1)-(3.3), 
with  q  denoting  the  velocity,  d  =  qx.  A  simple  manipulation,  (3.2)x2a+(3.3)x26 
yields  that  rj  :=  a2  +  b2  satisfies  the  mass  equation  (1.11), 

Vt  +  (vq)x  =  0,  rj:=a2  +  b2. 

Spatial  integration  of  (3.1)  implies  that  the  primitive  q  =  f  d  satisfies 

qt  +  qqx  =  vx,  Vxx  =  k{\a  +  ib\2  -  c)  =  k{i 7  -  c). 

We  conclude  that  ( r],q)T  is  a  smooth  solution  of  Euler-Poisson  system  (1.11)-(1.12), 
subject  to  initial  conditions  (r?o,<?o)T  =  (Po7Mo)T>  Thus,  by  uniqueness,  the  two  solu¬ 
tions  coincide,  (77  =  a2  +  b2,  q)T  =  (p,  u)T  for  t  <T*  and  hence  T*  <  T. 

Conversely,  assume  that  (p,u)T  is  the  smooth  solution  of  Euler-Poisson  system 
(1.11)-(1.12)  for  t  €  [0,T],  subject  to  initial  data  (po  >  0,zto)  such  that  po  =  +  ^o- 

By  assumption,  d  =  ux  is  smooth.  Consequently,  for  t  £  [0,T]  one  obtains  smooth 
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solutions  a  and  b  of  the  linear  transport  equations,  (3.2)  and  (3.3),  subject  to  initial 
conditions  (ao,5o)T-  Combining  the  two  equations  for  a  and  b  gives,  as  before 

(a2  +  b2)t  +  q(a2  +  b2)x  =  -(a2  +  b2)d. 

Subtracting  this  from  the  mass  equation  (1.11),  we  find  that  w  :=  p  —  (a2  +  b 2)  is  a 
solution  of  the  transport  equation 

wt  +  uwx  =  — wd , 

with  zero  initial  data,  w(x,  0)  =  0,  and  hence  w  =  0,  i.e.,  p  =  a2  +  b2 .  Differentiation 
of  (1.12)  then  yields 

(dxu)t  +  u(dxu)x  =  -k(p  -  c)  =  —k(\a  +  ib\2  -  c), 

which  shows  that  (dxu,  a,  b)T  is  a  smooth  solution  of  the  modified  system  (3.1)-(3.3) 
with  initial  data  (do,ao,b0)  for  t  G  [0,  T]  and  hence  T  <  T*.  □ 

Lemma  3.1  tells  us  that  the  regularity  of  the  extended  system  (3.1)-(3.4)  is  dic¬ 
tated  by  the  regularity  of  the  associated  ID  model  (1.11)-(1.12).  The  critical  threshold 
of  the  latter  was  analyzed  in  [ELT],  and  the  result  is  recorded  here  for  the  reader’s 
convenience. 

Theorem  3.1  (Global  regularity  with  sub-critical  data).  The  system  of  Euler- 
Poisson  equations  (1.11)-(1.12)  admits  a  global  smooth  solution  for  sub-critical  initial 
configurations  (po,uo)  with  JIR(p o  —  c)dx  =  0  such  that 

\dxu0(a)\  <  Jk(2p0(a)  -  c ),  Va  G  IR.  (3.5) 


The  precise  behavior  of  the  solution  of  (1.11)-(1.12)  is  given  along  the  particle 
paths 

(  -o  i  M(aj0)  .  (  /—+\  ,  /_oo(Po-c)d&  rr+\  to 

x(a,t)  :=  a  H - ^sin(vcu)-| - (1  —  cos (yckt).  (3.6) 

Vck  c 

If  (3.5)  holds  then  xa  ^  0,  and  we  have 

Do  , —  k  fx  (po  -  c)dx  . — 

p(x(a,t),t)  =  — u(x(a,t),t )  =  ug(a)  cos(v  ckt)-\ - 00  _ sin(vcfcf). 

xa{a,t)  V ck 

If  condition  (3.5)  fails,  then  the  solution  breaks  down  at  the  first  finite  time,  tc,  where 
xa(a,tc)  vanishes. 

Lemma  3.1  together  with  Theorem  3.1  yield 

Theorem  3.2  (Critical  threshold  for  modified  system).  The  limit  system  (3.1)- 
(3-4)  with  initial  data  (do,  do,  &o)T  with  J^R(a2  +  b2  —  c)dx  =  0  admits  a  globally  smooth 
solution  if  and  only  if 

|c?o(aOI  <  \J fc( 2ao(a)  +  26q(o)  —  c),  Va  G  IR.  (3.7) 

In  this  case  the  solution  of  (3.1)- (3.3)  is  given  by 


( a 2  +  b2)(x(a,  t),t) 


a2o(a)  +  bl(a) 
xa(a,  t) 


do  (a)  cos  (Vckt)  +  y/k/c(po(ce)  —  c)  sin  (\/ckt) 
xa(a,t) 
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4.  Local  convergence  rate.  According  to  Theorem  1.1,  the  modified  QEP 
solution  Ue  =  ( de,ae,be)T  remains  uniformly  bounded  together  with  couple  of  its 
derivatives,  independent  of  e.  Passing  to  the  limit  for  the  time  interval  t  <  Tq, 
the  QEP  solution  converges  to  a  limit  solution,  U,  which  is  governed  by  the  Euler- 
Poisson  equations  in  their  modified  formulation  (3.1)-(3.3).  We  want  to  estimate  the 
convergence  error.  To  prove  Theorem  1.2  let  us  consider  the  difference  We  :=  Ue  —  U 
which  is  governed  by  the  error  equation 

wte  +  qew :  =  (q-  qe)Ux  +  €-JW*xx  +  6-JUxx  +  Q{U\  t/e)  -  Q(U,  U). 

We  seek  an  error  bound  ||<92VF£||2.  A  straightforward  energy  method  yields 

ljt\\d2xwe\\2  =  (d2xdtwe,d2xw*) 

=  ( -  dltfdxw%dlw*)  +  ( dlu  -  q)dxu),d2xw*) 

+  \  (. Jd2x(Uxx  +  W*x)),d2xW*)  +  ( dl{Q{U £,  t/£)  -  Q(U,  17)),  32W£) 
=:  7i  +  /2  +  I3  +  /4. 


For  the  first  term  we  have  the  same  estimate  we  had  before  for  I\  in  (2.12) 


\h\  =  |(52(g£W£),d2W£)|  <  CW||17£||2||W£||2. 

We  now  turn  to  the  second  term 

h  =  (d2x((qe  -  q)dxU),  d2W*)  =  J2  (  l)  {dkx(qe  -  q)dl~kU ,  92W£) 

k—0  '  ' 


(4.1) 


k—0 


The  expressions  on  the  right  are  estimated  from  above  by  Sobolev  inequality 

|720|  =  |((9e  -  q)dlU,dlWe)\  <  ||<9|C7||oo -  <?||||W||2  <  C'onst.||17||4||We|||, 

\hi\  =  \2{{q^q)xdlU,dlW^)\ <  ||52t/||00||(g£-g),||||W£||2  <  Const. \\U\\3\\We\\l 
\h2\  =  \{dl{q*-q)Ux,d2xW*)\<  \\UX\\^U  -  qUW%2<C(mstm\2\\W%l. 


For  the  third  term  we  have 


\h\  =  e\(Jd2x(Uxx  +  W*x),d2xW*)\  =  e\(Jd2Uxx,d2xW*)\  <  e||t/||4||W£||2.  (4.2) 

Finally,  we  treat  the  forth  term,  /4.  We  first  note  that  difference  between  the  quadrat¬ 
ics,  Q(Ue,  Ue)  —  Q(U,  U)  must  be  of  the  form 

Q(Ue ,  Ue)  -  Q(U,  U)  =  L[U ,  ( Ue  -  17)]  +  R[(Ue  -  U),U%  (4.3) 

where  L  and  R  are  some  fixed  bilinear  form  (which  in  this  case  involve  the  constant 
k).  We  split,  accordingly,  /4  =  J4  +  iv4  and  we  find 

J4  :=  (02xL[U,  W%d2xW<)  =  E  (  fc  )  (L[dkxU,d2~kW%d2xW* 

k—0  '  ' 

and  we  find  with  the  usual  estimates 

2 

\M  <  Const.  E  ll^!7|U|W£||2_fc||W£||2  <  Const\\U\\3  ■  ||W£||2. 

k—0 
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A  more  careful  treatment  is  required  to  handle  I\4,  when  L  replaced  by  R, 
Ki  :=  ( dlR[W*,U%d2xW £) 

=  E  (  l  )  (R[dkxW%d2-kU%d2xW<)  =:J2R4k, 

k—0  '  '  k—0 

each  of  the  terms  on  the  right  is  upper  bounded  as  follows, 

\K40\  <  Con8t.\\We\\00\\Ueh\\We\\2  <  CoriSf.||C£||2||W£|||, 

\K41\  <  Const. Wd^W^WU^WW^h  <  Const. ||?7e||2||H/e||2, 
\K42\  <  C<m8t.\\We\\2\\Ue\\00\\We\\2  <  Const.\\U*\\2\\W%. 

Summarizing  the  upper  bound  we  have  so  far  we  conclude 


ft\\Weh  <  4U\ |4  +  C'1(||C/(-,t)||4  +  ||C£||2)  •  ||W£||2.  (4.4) 

According  to  Theorem  3.2,  the  particular  solution  U  and  its  derivatives  remain  uni¬ 
formly  bounded  in  time,  since  they  could  be  expressed  in  terms  of  dk(a2+b2),  dkd  and 
dku  which  are  all  periodic  in  time.  Let  C4  denote  the  upper  bound  || Z7*(- ,  i) ||4  <  C4. 
By  Gronwall’s  inequality  the  error  does  not  exceed 

||t/£(-,f)-C(-,f)||2  (4.5) 

<  ||t/e(-,0)  -u(;  0)||2  +  C4e^  0exp{JT  Cl(C4  +  ||t/£(-,  r||2)dr  |  ds. 

We  recall  that  for  the  short  time  interval,  t  <  To,  we  have  by  Theorem  1.1, 
||tC(-,f)||2  <  2 1 1 Z7o 1 1 2  and  (4.5)  implies  convergence  rate  of  order  0(e)  asserted  in 
Theorem  1.2. 


5.  Global  convergence  of  semi-classical  limit.  The  error  equation  (4.4) 
shows  that  during  the  local  existence  interval,  [0,To],  the  QEP  solution  cannot  grow 
too  fast  as  predicted  by  the  pessimistic  Ricatti  bound  (2.13), 

|l|c£||2<c0||t/£||2. 

Indeed,  Ue  must  stay  within  the  0(e)  neighborhood  of  the  limit  solution  [/(•),  and  in 
particular, 


\\U4;t)h  <  ||[/(•,f)||2  +  ||^/£(•,f)-[/(•,f)||2  <  C2  +  C4e  f  e/T‘=s  CRCi+\\U‘(;rh)dr  dfj 

J  s- 0 


We  want  to  extend  this  argument  beyond  the  local  existence  interval  [0,7b].  The  error 
bound  (4.4)  is  valid  for  all  f’s,  linking  the  growth  of  the  error  e(t)  :=  ||C£(-,  <)  —  [/(•,  f)||2 
to  the  growth  of  \\U(-,  i)||2, 


—  e(t)  <  C4e  +  Ci (C4  +  ||C£||2)  •  e(t). 


Together  with  the  upper  bound  ||£/£||2  <  \\U\\2  +  \\Ue(-,t)  —  U(-,t) ||2  <  C2  +  e(t)  we 
conclude 

j2e{t)  <  C4e  +  Cl(C4  +  C2  +  e(f))  •  e(t)  <  C[e  +  e(t)  +  e2(f)],  e(0)  =  0. 
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This  gives  the  inequality  for  Z(t)  =  e(t)e  ct, 


jf  Z(t)  <  Cee~ct  +  CeCtZ2,  Z{ 0)  =  0. 

Multiplying  the  above  inequality  by  (1+gZy  with  6  chosen  as 

f  eCT  ) 

9  =  max  |4(eCT  -  1),  A^CT  _  ^  j , 

we  obtain 


6Z' 

(i  +  ezy 


<  C9ee~ct  +  CO^e0*. 


Integration  over  [0,  t]  gives 


1 

1  +  9Z 


>l-e0(l  -  e~Gt)  -  6~l{ect 


1)> 


1 

2 


for  e  <  £q  with  e0  =  16(ecT_1) 2  and  0  chosen  above.  Therefore  we  have  Z  <  1/0,  i.e., 


e(t)  <  min 


oCT 


A(eCT  -  1)  ’ 


4e(. 


CT 


This  shows  that  for  any  T  >  0,  there  exists  a  sufficiently  small  £0  ~  e  CT,  such  that 
for  e  <  e0  we  have  \\Ue(-,t)  —  U(-,t) H2  <  Constre  and  Theorem  1.3  follows. 

6.  Validity  of  the  WKB  expansion.  According  to  Theorem  1.3,  the  NLSP 
equation  (1.8)  admits  a  global  solution  with  modified  variables,  U£  =  ae ,  be)  such 
that  Ue  =  U  +  0(e)  where  U  satisfies  the  limiting  modified  system  (2.11).  The  latter 
is  an  hyperbolic  system  so  that  U(-,t)  is  in  Hs  if  [/(•,  0)  is.  We  want  to  show  that 
the  NLSP  solution  ip£  =  A£exp{ii f»£/e},  realized  in  terms  of  the  3-vector  of  modified 
variables  U£  :=  [4)^,  a£,  b£],  admits  the  following  asymptotic  expansion 


N 

U*  =  Y,Uj£j  +o{eN), 


j=o 


Uj  g  L°°([0, T],  (IR)). 


•  Leading  order 

From  Theorem  1.3  we  know  that  U£  are  bounded  in  L°°([0,  T\,  1L2(IR)).  Using 
the  governing  equation  (2.11),  [//  are  bounded  in  L°°(0,  T;  L2(1R));  in  fact, 
according  to  (2.11),  Ue  is  given  as  the  sum  of  three  terms:  an  H 1  x  H 1  term 
q£U£,  a  vanishing  L 2  term,  ~  tJU£x  and  an  H 2  quadratic  term,  Q(U£,U£), 
which  are  compact  in  L2.  It  follows  that  a  subsequence  of  U£  converge  uni¬ 
formly  in  L°°([0,  T],  L2(IR))  to  Uq,  which  are  solutions  of 

dtUo  +  qodxUo  =  Q(Uo,Uq),  Uq  =  [do,  ao,ba],  do  =  dxqo  =  dxx$o 


subject  to  Hs  initial  data 


Uq{-,  0)  =  lim  U£(-,  0)  G  Hs. 

e— >0 

Hyperbolicity  implies  that  Uo(-,t)  remains  in  Hs  and  the  uniqueness  of  limit 
solution  for  this  system  ensures  the  convergence  of  the  whole  sequence,  U£  — > 

U0. 
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First  order 

To  justify  the  first  order  term  in  the  expansion,  we  consider  the  iP_2-initial 
data 


Hi(-,0)  =  lim 

e — *0 


He(-,0)-H0(-,0) 

e 


Setting  We  :=  ( Ue  —  Ho)  A,  and  carrying  out  the  energy  estimates  along  the 
lines  of  §4  leads  to  an  upper  bound 


ll^eM)||s-2<.F(||t/oy,  t  <t, 


with  a  nonlinear  dependence  IF,  depending  on  {||Ho(-,f)||s}o<t<T-  Hence 
We  is  bounded  in  L°°([0,  T],  iP-2),  and  Wte  is  bounded  in  L°°([0,  T],  FP-4). 
As  before,  a  subsequence,  We  converges  strongly  in  L°°([0,  T\,  iP-4)  to  a 
limit  function  U\  =  [q±,  ai,bi\.  This  limit  solves  the  following  first  variation 
problem 


1 

to 

a- 

o 

2kao 

2kbo 

dtUi  +  qodxUi  +  qidxU0  —  Q(Uq,  Uq)  + 

o  do 

-fso 

2^0 

0 

Hi,  (6.1) 

0 

—  2 

with  initial  data 


Z7i(0)  =  lim 

e — *0 


ue( o)  -  H0(Q) 

e 


The  hyperbolicity  of  (6.1)  implies  that  U\(-,t)  remains  in  iP~2  and  compact 
in  iP-4.  The  uniqueness  of  solution  to  this  problem  ensures  the  convergence 
of  the  whole  sequence  We  to  U\  £  IP-4. 

Higher  order  terms 

Assume  that  we  have  already  an  asymptotic  expansion  to  the  order  N 


N 

C/e  =  ^C/,P+o(eJV), 


3=0 


with  Uj  £  L°°([0,  T\,  Hs  2j).  Let  We  =  W ^  denote  the  error 


We 


Ue  -  U€ 

eN+ 1  ’ 


3=0 


Inserted  into  (2.11),  we  find  that  We  satisfies 

dtW*  +  qfdxWf-  +  {q'  -  cf  )dxUe  =  e-JWlx  +  ^  JUexx  +  tP)  (6.2) 

where  iP  =  H eN  is  given  as  the  sum  of  bilinear  form,  consult  (4.3) 

HeN  =  L(U\  We)  +  R(Ue,  We )  +  eN+1Q(We,  We). 

One  can  show  that  HeN{-,t)  is  bounded  ||-ff/v||s-2jv-2  <  Const. || We ||^_2Ar_2. 
Assume  that  the  initial  data  satisfy 


||!Fe(0)||s_2AT-2  <  Const.., 
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then  (6.2)  implies  that  the  scaled  error,  e(t)  :=  ||W(-, f||s_2Ar-2>  satisfies 
e(t)  <  C2s  +  C\e(t)  +  C2e2(t),  which  in  turn  yields 

\\W^;t)\\s-2N-2  <  Const,  t  <T  ~  N  In 

e 

We  conclude  that  \\Ue  —  Uje^\\a-2N-2  <  0(eN+1). 

•  WKB  expansion 

The  above  results  show  that  we  have  obtained  the  formal  expansion  of  Ae 
and  4>e  to  an  arbitrarily  high  order.  To  recover  the  usual  WKB  expansion,  it 
suffices  to  equate  the  e-powers  of  the  same  order  in  the  identity 

oo  _  /  oo  \  /  oo  \ 

E  =  E  exP  *  E  • 

j— 0  \ra=0  /  \  m— 0  / 
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